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The homoepitaxial growth of initially flat surfaces has so 
far always led to surfaces which become rougher and rougher 
as the number of layers increases: even in systems exhibiting 
"layer by layer" growth the registry of the layers is gradu- 
ally lost. We propose that pulsed glancing-angle sputtering, 
once per monolayer, can in principle lead to layer-by-layer 
growth that continues indefinitely, if one additional parame- 
ter is controlled. We illustrate our suggestion with a fairly 
realistic simulation of the growth of a Pt (111) surface, cou- 
pled with a simplified model for the sputtering process. 



When depositing atoms on a surface, one often wants 
control of the growth morphology. At the most primi- 
tive level, when depositing atom X on a low-index, flat 
surface of the X crystal, it would be nice to be able to 
ensure that the resulting surface remained flat! 

This simple goal has often been difficult to achieve in 
practice, especially in metals. At low temperatures, one 
observes three-dimensional growth: inter-layer mobility 
is low, and the second layer starts growing as soon as 
the flrst layer gains any substantial coverage. At inter- 
mediate temperatures, one observes what is called two- 
dimensional or layer-by-layer growth. Layer-by-layer is 
used to describe systems with oscillations in some mea- 
sured property (antiphase scattering of He [Q , RHEED 
1^, low-energy electrons ||], or X-rays Q), decaying 
slowly as the number of layers increases. For a few ma- 
terials under special conditions, this decay can be quite 
slow Q, but under most circumstances it decays over a 
few tens of layers. At higher temperatures, for slightly 
mis-cut surfaces, one can have a step-flow regime which 
typically exhibits quite stable layered growth. 

In this paper, we consider the question of how one 
might achieve indefinite layer-by-layer growth: a mode- 
locked state |§] where the surface irregularities due to the 
growth process would remain bounded and oscillations 
in the properties would continue forever. Such long- 
range order despite the noise (random fluctuations in 
deposition, nucleation, and growth) would be separated 
from the traditional decaying layer-by-layer growth by 
a phase transition. The key is to periodically force the 
system, pulsing in synchrony with the deposition of each 
monolayer (keeping in phase by using one of the real- 
time oscillatory measurements described above). Pulsing 
the temperature, pulsing the deposition rate, and pulsing 
with an ion-assisted anneal have been used to good effect 
experimentally. |jl| We argue that pulsing these quanti- 
ties cannot lead to indefinite flat growth, but that pulsed 
sputtering can yield layer-by-layer growth forever. 



The dashed line in figure n shows a typical ther- 
mal growth on a surface. It is a numerical simu- 
lation of Pt/Pt(lll) at 130 K, grown at one mono- 
layer/second, with parameters determined using Effec- 
tive Medium Theory and available experimental infor- 
mation as described in B. The figure shows / — 
(^^q(— l)*(0i+i — 6*^)) as a function of time, where 
6i is the fractional coverage in the ith layer; this mea- 
sure corresponds closely to what is measured in antiphase 
scattering probes used in the experimental systems. No- 
tice that the signal vanishes whenever the number of 
atoms in odd and even layers become equal. Notice that 
the peaks decay as more layers are deposited: as the sur- 
face begins to span several layers, the surface morpholo- 
gies at integer and non-integer monolayer coverages be- 
come indistinguishable. One must note that, for the pa- 
rameters we simulate, the surface does not grow wildly 
rough even for thermal deposition, where we observe four 
oscillations in /. 
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FIG. 1. Layer by layer growth, with and without pulsed 
sputtering. The main plot shows the simulated antiphase 
intensity 7; the oscillations with a traditional thermal de- 
position (dashed line) decay rapidly. Pulsed sputtering of 
A — 25% of a monolayer is done after each deposition of 1.25 
monolayers, yielding sharp jumps in the curves. (Here the 
sputtering phase ^ = 0.20.) The lower curve in the inset 
shows the peaks of the jumps for the first hundred monolay- 
ers deposited: we expect the layers to stay flat indeflnitely, for 
the correct choice of /i. The upper inset curve is the signal 
from an initially rough sample prepared by first depositing 
several layers without sputtering. When the sputter-deposit 
sequence is started, we observe that the surface regains a fiat 
interface. 

This decaying layered growth is usually described by 
theories focusing on the nucleation and growth of islands 
on the surface. ^ Most of the pulsed attempts to improve 
the stability of layered growth [y have been motivated 
by the nucleation and growth theories, and have deliber- 
ately increased the nucleation rate at the onset of a new 
layer (where depositing on top of existing islands is not 
a concern), while reducing it thereafter. 



There is another school of theoretical models, which 
focus not on individual islands but rather on the effect 
of fluctuations in the deposition rate and the role of dif- 
fusion within a continuum description for the height of 
the surface. These models predict that the random fluc- 
tuations in deposition will always overwhelm the avail- 
able diffusive mechanisms for retaining a flat interface, 
on sufficiently long length and time scales. |lQ] There 
are several mechanisms and models [0| for this diffu- 
sive smoothing of surfaces, but they all predict that the 
surfaces will eventually become rough: the noise is inde- 
pendent of wavelength, and the diffusion becomes feeble 
at long wavelengths. However unlikely it is to nucleate 
on top of an existing island, diffusion cannot transport 
the extra atoms from one region of the surface to an- 
other fast enough: eventually the extra atoms in one re- 
gion will nucleate extra layers. Pulsing the temperature, 
pulsing the deposition rate, or pulsed annealing with an 
ion beam only changes the effective diffusion rates on the 
surface, and does not fundamentally alter this conclusion: 
we need a non-diffusive mechanism. 

How can we smooth the surface in a more effective 
way? High-energy atomic beams incident at glancing 
angles to the surface are a known way of generating flat 
surfaces |13): the beam preferentially sputters atoms off 
the mountains and hills. Especially for groups using en- 
ergetic beams for growing surfaces [Q , it would seem nat- 
ural to try to use a pulsed sputtering mechanism. Con- 
sider starting with a flat surface with an initial deposition 
of f + /i monolayers. We then repeatedly sputter off A 
and deposit 1 -I- A monolayers, so as to always start the 
sputtering at an integer plus jjl monolayers coverage. If 
a surplus of atoms is deposited onto a region of the flat 
surface, there will be a surplus of adatoms at the time 
of sputtering, and thus the sputtering will remove ex- 
tra mass from the region. This non-diffusive mechanism, 
being independent of wavelength, should in principle be 
able to produce indefinite layered growth. Q 

We decided to test these ideas with a simple model. 
Building on our well-characterized, g physically realistic 
solid-on-solid kinetic Monte-Carlo model for thermal de- 
position and growth on Pt (1 1 1), we implemented a prim- 
itive, simplistic model for glancing-angle rotating-beam 
sputtering. |I3] We wanted atoms in an intact layer to 
be immune from sputtering, solitary atoms and edges of 
small islands to be sputtered at high rates, and pit edges 
to be at least partially shielded. We sputter equally from 
each of the six directions lying along atomic rows. An 
atom is immune from sputtering from a given direction 
if it is shadowed by another atom upwind in the same 
row and the same monolayer, within a distance L — b 
interatomic distances. Figure ^ shows that this model 
can produce smooth growth, seemingly forever. 

This sputtering model we imagine might correspond 
roughly to a beam at an angle arctan(I/L) = 11°. Our 
model is likely over-optimistic in that intact monolay- 
ers and shadowed atoms are protected completely. It is 
using a pessimistic angle of attack: angles from 1° to 



15° have been used for smoothing rough surfaces [ [l3[ , 
and a more glancing angle in our simulation would pro- 
duce much smoother surfaces. This kind of simplistic 
model is particularly useful in studies of the qualitative 
features and feasibility of a new method, and indeed it 
immediately uncovered an important issue neglected in 
our discussions so far. 
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FIG. 2. Critical pit. This is a snapshot of our simulation 
with A — 0.25, and /i = 0.1, after 41 layers have been de- 
posited. Lower layers are darker grey (pits). The inset shows 
a plot of the percentage of the system size occupied by the 
largest island and largest pit, as a function of how many mono- 
layers have been deposited. The arrow in the inset shows the 
time of the snapshot. Until 39 layers, all we see are small pits 
and islands which appear and disappear. At this point, we see 
a large pit form (circled), which grows irregularly thereafter. 
The oscillations in the anti-phase intensity (as in figure hi) die 
out slowly after the pit nucleates. 

Figure shows a surface grown with our sputtering 
schedule with a different value of fj,, after growing 41 
monolayers. Notice the pit. Up until this point in the 
simulation, we saw excellent layer-by-layer growth, simi- 
lar to that shown in figure 0; subsequent to this frame the 
oscillations die away. We interpret this behavior as the 
nucleation of a critical pit, analogous to critical droplets 
at first-order phase transitions. |14[| 

Consider what happens to a pit of radius i?„ under a 
cycle of depositing 1 + A monolayers and sputtering off A. 
During the deposition, atoms falling on the upper layer 
will typically nucleate new islands, which grow and merge 
to raise the height by one. Atoms landing inside the pit 
will stick to its outer edges; the pit will fill in (given 
a large Ehrlich-Schwoebel barrier) at about the time a 
whole monolayer is deposited. However, the region of 
the former pit will remain depressed, since new islands 
start nucleating and growing only after it fills in. As a 
zeroth approximation, a new pit one level higher of the 
same radius will exist after one cycle of deposition and 
sputtering: Rn+i ^ Rn- 

Consider a fiat step on the surface — interpretable 
either as a pit or an island of infinite radius. Under one 
cycle of deposition and sputtering, there is no reason to 
expect that the attachment at the step edge will balance 



the sputtering. (That is, our model has no symmetry 
between pits and islands.) One expects the step edge 
to move by a distance A, where we define positive A to 
represent the growth of the pit (lower terrace). To a first 
approximation for large radius, we expect the pit area 
after one cycle on average to change: i?„+i = Rn + A. 

Now, since the edges of small pits are partially 
shielded, and islands are more exposed to the sputtering, 
the net effect after one cycle is to remove less material in 
existing pits, and more near existing islands. The smaller 
the pit radius Rn, the more protected is the pit, and the 
higher is the net deposition after the entire cycle. To a 
second approximation, we expect that the pit will be re- 
formed at a new radius Rn+i = Rn + A — E/i?„. The 
term — S/i?„ represents the physics of the self-shielding 
for small pits which tends to make them shrink; it also 
makes small islands shrink (negative radius). This is the 
term found for detachment-limited coarsening for islands 
on surfaces [^ , is the first term in a Taylor series in the 
curvature of the edge, and can be derived with a simple 
geometrical argument based on shielded sites. 
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FIG. 3. Large pits grow, small pits decay, for A — 0.25, and 
fj, — 0.15. Shown is the evolution of the area Ap of a pit for 6 
different initial pit sizes A, each 10% of the initial area of the 
simulation, averaged over 100 runs each (except for size 100, 
with 85 runs). We follow the evolution of the pit by measuring 
its size every time there is an integer number of monolayers 
down. The thin curves are equation hi fit to the data (if a 
is the lattice spacing, then A ~ 0.09a, and E ~ 0.9a^, with 
an initial shrinkage of the island radius of about 6%.) The 
critical pit size Re is estimated to be 280 ±60 in area, or about 
ten lattice constants in radius. 

Finally, there are the stochastic fluctuation in the pit 
radius. We expected that the dominant source of fluctu- 
ations would be the fluctuations in the number of atoms 
deposited within the area of the pit, which should scale as 
the square root of its area: hence the radius fluctuations 
will be independent of radius. 

Rn+l=Rn + A-j:/Rn + n^n- (1) 

ri £,n in equation m is the noise term: ^„ is a random vari- 
able with mean zero and standard deviation one, and CI 
gives the strength of the noise. Measuring these fluctua- 
tions directly, in the interesting range < /i < 0.2, we've 



verified that they are indeed roughly independent of R 
and of width Q ~ Q.65a. Figure |3| shows the results of 
a fit of the average shrinking and growing of pits, using 
equation n^ without noise (fJ = 0). 

The layer-by-layer growth in our model ends when the 
fluctuating noise produces a pit of the critical size 



(i?e) - S/A, 



(2) 



after which the pit grows by itself to macroscopic size. 
Our equation m can be thought of as a thermal randonr- 
walk in radius, with step size 17, temperature T — 2f2^, 
and potential S log(i?) — (i?— 1)A: |Q the critical radius 
is the local maximum Vmax in the potential. One can 
solve a continuum approximation to equation |l| for the 
rate of formation of large pits, per density of pits of size 
R^l: 



J = A(2A/f72)2E/o^/p[^ _^ 2S/f72, 2A/f72] 
= A {(2E/(ef72))2SA2'/r[l + 2T./n'\ 2A/f}2]| , 
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where F is the incomplete gamma function. 

The last of equations g shows the connection with tra- 
ditional critical droplet theory. Q,|5[ Here the term in 
curly brackets is a prefactor, A is a bound on the velocity 
at which one could cross the barrier, and e"^™"^/-^ is the 
Boltzmann probability of sitting at the critical radius. 

How can we grow layer-by-layer forever? Clearly, we 
wish to set A to zero, imposing a long-wavelength sym- 
metry between islands and pits sending the nucleation 
rate J to zero. All three constants A, S, and f7 in equa- 
tion (Q) will depend on temperature, deposition rate, 
sputtering angle, other adsorbates on the surface, the 
fraction A sputtered, and the point /i during the depo- 
sition of a monolayer that the sputtering occurs. If by 
varying any of these parameters we can set A = with- 
out making E < 0, we ought to suppress the nucleation 
altogether, and sustain layered growth indefinitely. 

Fitting to simulations like those shown above in fig- 
ure 0, we have measured the critical pit size R^ A, and 
S as functions of ^. We find A, E, and an initial is- 
land shrinkage all fit well to the form asin(27r/i -\- <j)) +b. 
Direct measurements of the critical pit size show a diver- 
gence where our sinusoidal interpolation for A changes 
sign (figure ^). [0 Above /Xc, where A < 0, we expect 
all pits to be stable and large islands to be unstable, flq ] 

Can we show that the surfaces remain flat near ^c^. 
The sputtered simulation shown in figure |l| was done at 
/i — 0.20 (near jjLc — 0.24 where the lifetime of fiat growth 
diverges); the inset shows that the oscillations are per- 
sisting as long as we have simulated. Equation (0) pre- 
dicts the rate of formation of large pits J{^ = 0.20) to 
be 1.4 X lO""*^ times the density of pits of size one. For 
values of /i far from /ic the oscillations decay rapidly: at 
/i = 0.7 the oscillations die roughly as they do for a ther- 
mal growth without sputtering (although the rms rough- 
ness for the sputtered surface is much smaller than that 
for a thernrally grown, unsputtered surface even away 



(3) 



from jic)- Figure ^ also shows our prediction of the nu- 
cleation rate J for large pits: proportional to the inverse 
of the lifetime for flat growth. 
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FIG. 4. Critical pit size Re (in units of the lattice constant 
a) and large pit formation rate J (in units of nucleated large 
pits per monolayer per vacancy) as a function of sputtering 
phase /I. The points with error bars are eyeball estimates from 
plots like figure H. The upper curves are from equation eI with 
sinusoidal fits to A(/i) and E(/^). Because A changes sign 
at fic ~ 0.24, the critical pit size diverges there; for fi < fic 
large pits are unstable (solid curve), and for fi > fic large 
islands become unstable (dashed curve). When the critical 
pit size diverges, the rate of large pit formation J goes to 
zero (dot-dashed curve). 

Finally, in perhaps the most convincing demonstration 
that our method is working, note back in figure H the sec- 
ond plot in the inset. It shows a rough surface becoming 
flat once we start our pulsed sputtering schedule! 

In trying to grow flat layers, why not go for perfec- 
tion? We argue that the methods used heretofore to grow 
flat surfaces (pulsing temperature, deposition rate, ion- 
assisted diffusion, ...) cannot compete in the end with 
the stochastic noise in the deposition rate. We claim that 
pulsed sputtering, smoothing once per deposited mono- 
layer, can in principle yield layer-by-layer growth oscil- 
lations that last forever, provided that one parameter is 
tuned to a critical value. 
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